wo 2004/077359 



PCT/SG2003/000043 



26 



Claims 

1. A method for extracting cerebral ventricular system information from images of one 
or more cerebral ventricular regions, the cerebral ventricular system comprising a 
third ventricle (V3), a fourth ventricle (V4), a left lateral ventricle (VLL) having a 
body (VLL-B) and an inferior (temporal) horn (VLL-I), a right lateral ventricle (VLR) 
having a body (VLR-B) and an inferior (temporal) horn (VLR-I), an anterior 
commissure (AC), a posterior commissure (PC), and a midsagittal plane (MSP), the 
method comprising the steps of: 

1) defining multiple regions of interest (ROI) in -the images; 

2) defining seed points within each ROI; 

3) growing images of ventricular regions while correcting for leakages into 
extraventricular space; and 

4) connecting the ventricular regions grown. 

2. A method according to claim 1, \<4ierein the steps of defining multiple ROIs, 
defining seed points and growing images are applied first to one or more images 
of a third ventricle (V3), then to one or more images of the fourth ventricle (V4), 
and then to the left and right ventricles (VLL and VLR) for controlling leakage 
and connections. 

3. A method according to claim 1 or 2, wherein the method is performed from one or 
more medical imaging modalities. 

4. A method according to any one of the preceding claims, wherein the step of 
defining multiple ROIs comprises defining each ROI in a predetermined plane. 

5. A method according to claim 4, whereia the step of defining multiple ROIs 
comprises selecting said predetermined plane to be on a coronal orientation to 
constitute a coronal slice. 

6. A method according to claim 5, wherein the step of selecting said predetermined 
plane to be on a coronal orientation for assessment of an image of the body of the 
left lateral ventricle and/or the body of the right lateral ventricle comprises the 
steps: 

(1) selecting the coronal plane (VAC) passing through the anterior 
commissvire(AC); 

(2) selecting the initial rectangular ROI on the VAC laterally between the co- 
ordinates [MSP, MSP+nl] for the body of the left lateral ventricle, and [MSP- 
nl, MSP] for the body of the right lateral ventricle, and dorsally between 
[AC, AC+n2], where nl and n2 are constants; 

(3) altering the ROIs for subsequait processing based on desired histogram 
distribution within the ROIs. 

7. A method according to claim 6, wherein each aspect ratio of the ROIs is constant. 
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8. A method according to claim 7, wherein the sizes of the ROIs are reduced 
laterally and dorsally. 

9. A method according to any one of claims 6 to 8, wherein nl is aroimd 25nmi. 

10. A method according to any one of claim 6 to 9, wherein n2 is aromid 35mm. 

11. A method according to claim 5, wherein the images of the one or more cerebral 
ventricular regions comprises one or more coronal slices comprised of pixels, and 
wherein the step of selecting said predetermined plane to be on a coronal 
orientation for assessment of an image of the inferior hom of the left lateral 
ventricle (VLL-I) and/or the inferior hom of the right lateral ventricle (VLR-I) 
comprises the steps: 

1) selecting coronal slices between the coronal plane passing through the anterior 
commissure (VAC) and the coronal plane passing through the posterior 
commissxire (VPC); 

2) drawing a horizontal profile at the level of the inferior pixel of the image of the 
third ventricle(V3); 

3) identifying the temporal bone signal in the horizontal profile as the lateral-most 
high signal due to bone marrow fat; 

4) measuring the distance firom the temporal bone to the midsagittal plane (MSP); 

5) dividing the distance in half to define a midpoint and considering this midpoint 
as point "a"; 

6) drawing an ROI having the dimensions n6 x n6 in reference to the midpoint 
"a" so that the coordinates are (a, a-n7 (a, a+n7) laterally and (a, a+n6) 
inferiorly, where n6 and n7 are constants; and 

7) maintaiiiing the ROIs for subsequent processing unchanged. 

12. A method according to claim 1 1, wherein n6 is around 20mm. 

13. A method According to claim 1 1 or 12, wherein 7 is around 10 mm. 

14. A method according to claim 1 wherein said ROI for the third ventricle (V3) is set 
on a sagittal plane by: 

1) selecting the midsagittal plane (MSP); 

2) setting an initial rectangular ROI antero-posteriorly between the co-ordinates 
[AC, PC] and dorsally between [AC - n8, AC + n9], where n8 and n9 ai« 
constants; and 

3) altering the ROIs for subsequent processing based on desired histogram 
distribution within the ROIs. 

15. A method according to claim 14, wherein the sizes of the ROIs are reduced by 
moving the superior boimdary of the ROI inferiorily. 



16. A method according to claim 14 or 15, wherem n8 is aroimd 10mm. 
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17. A method according to any one of claims 14 to 16, wherein n9 is around 35mm. 

18. A method according to claim 1 where said ROI for the fourth ventricle (V4) is set 
on sagittal plane by: 

1) selecting the midsagittal plane (MSP); 

2) setting a rectangular ROI posteriorly between [PC, PC + nlO] and dorsally 
between [PC, PC - nl 1], where nlO and nl 1 are constants; 

3) limiting the superoposterior part of the ROI by drawing a line startmg at nl2 
posterior to PC to the posteroinferior comer of the ROI to bisect the ROI 
obliquely, where nl2 is a constant; 

4) limiting the inferoposterior part of the ROI by drawing the line from the 
inferoanterior comer of the ROI to the midpoint of the posterior limit of the 
ROI; and 

5) altering the ROIs for subsequent processing based on desired histogram 
distribution within the ROIs. 

19. A method according to claim 18, wherein the sizes of the ROIs are reduced 
posteriorly, 

20. A method according to claim 18 or 19, wherein nlO is around SOmm. 

21 . A method according to any one of claims 1 8 to 20, wherein nl 1 is around 45mm. 

22. A method according to any one of claims 18 to 21, wherein nl2 is around 5 mm. 

23. A method according to any one of claims 6, 14, or 18, wherein the desured 
histogram distribution within said ROIs the peaks corresponding to white matter 
(WM), grey matter (GM), and cerebrospinal fluid (CSF) in said histogram are 
identifiable with any existing peak detection algorithm. 

24. A method according to any preceding claim, wherein the step of defining seed 
points for within an ROI for VLL-B and VLR-B comprises the steps: 

1) sampling the ROI horizontally, starting from AC+nl3 every nl4 distance, 
where nl3 and nl4 are constants to form ample line segments; 

2) calculating a profile along each sample line segment; 

3) determining the longest CSF segment from said profile; and 

4) placing the seed point in the middle of the segment. 

25. A method according to claim 24, wherein nl3 is aroimd 10mm. 

26. A method according to claims 24 or 25, wherein nl4 is around 5 mm. 

27. A method according to claim 11, wherein the step of defining seed points for 
VLL-I and VLR-I comprises the following steps: 
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1) setting the ROIs for the VLL-I and VLR-I on nmltiple coronal planes between 
VAC and VPC, wherein the planes are selected and processed starting fixMn 

2) locating the CSF region on the coronal slice between VAC and VPC having 
the largest CSF region; and 

3) placing the seed point at the geometric centre of said CSF region. 

28. A method according to any preceding claim, wherein the step of defining the seed 
point for V3 comprises the steps: 

1) calculating a profile along the line between the anterior and posterior 
commissures (AC-PC); and 

2) determining a pixel on the AC-PC line, whose gray level is closest to the CSF 
mean value to constitute the seed point. 

29. A method according to claim 18, wherein the step of defining the seed point for 
V4 comprises the steps: 

1) sampUng the ROI horizontally, starting fi-om PC-nl5 every nl6 distance, 
where nl5 and nl6 are constants to generate a sample line segment; 

2) calculating a profile along each sample line segment; 

3) calculating the length of the CSF in each profile and comparing the calculated 
length to the previous one; 

4) when the calculated length decreases for at least nl7 subsequent slices, select 
the middle of the longest CSF segment as the seed point, where nl7 is any 
constant. 

30. A method according to claim 29, wherein nl5 is around 10mm. 

3 1 . A method according to claims 29 or 30, wherein nl6 is around 2 mm. 

32. A method according to any one of claims 29 to 31, wherein nl7 is around 2mm. 

33. A method according to any preceding claim, wherein each lateral ventricle is 
grown in 3D space on coronal slices, slice by slice. 

34. A method according to any preceding claim, wherein the step of growing the 
image of the third ventricle V3 comprises the steps: 

1) subdividing V3 into four subregions 1, 2, 3 and 4 by the planes passing VAC, 
VPC, AC-PC, subregion 1 containing the most anterior part of V3 and being 
separated firom subregions 2, 3 and 4 by the VAC and the AC-PC planes, 
subregion 2 containing the most dorsal part of V3 and being separated fi-om 
the subregions 1, 3 and 4 by the AC-PC and the VPC planes, subregion 3 
containing the most posterior part of V3 and being separated &om subregions 
1, 2 and 4 by the AC-PC and the VPC planes, and subregion 4 containing the 
most ventral part of V3 and being sq)arated firom the subregions 1, 2, and 3 by 
the AC-PC and theVAC planes; and 
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2) growing V3 in three dimensions, wherein subreeion i ic * • , 

^.^'c^^pSS fj^^''- in eS;t«Lt 

belo^ the PC in toSoTSS^ ^ """""^ « '-Sion 

41 . A method according to claim 40, wherein n 1 8 is around 4 



mm. 
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42. A method according to claim 35, wherein during the growth of the V4 region 
the step of correcting for leakages dorsopostetiorly tibrough the superior medullary 
velum to the cistema ambiens, comprises: 

1) hindering the number of foreground pixels of the current axial sUce from 
increasing by 50% more than that in the previous slice in the inferior direction; 
and 

2) maintaining the distance between the centre of gravity of the foreground 
region of V4 and the MSP less than around 2 mm 



43. A method according to claim 35, wherein during the growth of the V4 region 
the step of correcting for leakages ventroposteriorly through the inferior medullary 
velum to the cistema magna, comprises maintaining the width of the foreground 
region of V4 at the lateral recesses less than around 2 mm. 

44. A method according to any preceding claim, further comprising repeating the 
region growing with narrower growing ranges after detection of leakage. 

45. A method according to claim 44, wherein the step of repeating the region 
growing with narrower growing ranges comprises narrowing the intensity range 
and/or the spatial range of the growing region. 

46. A method according to claim 45, wherein the narrowing of the intensity range 
comprises decreasing the value of the scaUng factor s, in the equation defining the 
growing range of intensity [ -s*( -m), +s*(M- )] , wherein m is the 
nmumum, M the maximum and the mean values of the CSF range and s is a 
scaling factor between 0 and 1. 

47. A method according to claim 46 wherein the narrowing of the intensity range 
is performed iteratively by decreasing the value of s, followed by performing the 
region growing again and checking whether the leakage is avoided. 

48. A method according to claim 1, wherein the step of connecting the ventricular 
regions VLL-B and VLL-I comprises: 

1) locating the leftmost inferior pixel of [VLR-B, Pi ] where Pi is the left most 
inferior pixel of VLR-B; 

2) growing VLR-I on a coronal slice in a posterior direction; 

3) checking whether the number of foreground region is equal to or less tiian 1; if 
not, go to 4); if yes, 

• find the minimum rectangular area containing all the foreground pixels 
of the previous coronal slice in anterior direction; 

• pass the rectangular area to the current coronal slice C : 

CUT * 
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• find a seed point for the grey matter at the rectangular area, and extract 
the grey matter through region growing; if the seed point cannot be 
found or the grey matter is absent, go to 7) for brute-force connection; 

• trace the leftmost supraior boundary of the grey matter; 

• add the traced pixels to the grown region, and go to 5); 

4) checking whether the number of foreground region(s) is more than 1 on ; 
if not, go to 5); if yes, comiect the foreground regions by the following steps: 

o locate the P2 and P3 pixels on the foreground regions, where P2 is the 
leftmost inferior pixel of the upper foreground region and P3 is the 
leftmost superior pixel of the lower foreground region; 

o trace the boxmdary of the grey matter between P2 and P3; 

o add the traced pixels to the grown region; 

5) checking whether the current grown region is connected with the grown region 
on the previous coronal slice in anterior direction; if not, connect them; 

6) checking whether Pi is on the current coronal slice; if not, go to 2); if yes, 
terminate the current connection process; 

7) enforcing the brute-force connection when either fhs seed point is not found or 
the grey matter is absent; and 

8) forming a cone from Pi and the last grown region to complete the connection. 

49. A method according to claim 1 , wherein the step of connecting the ventricular 
regions V3 and V4 after finishing segmentation of V3 and the subregion 1 of V4 
comprises: 

1) extracting the connection in the superior direction by: 

(a) locating the superior connection starting pixel (Qi) at the posterior 
inferior part of V3; 

(b) locating the most supraior pixel of subregion 1 of V4 (Q2) 
assuming the location is inferior to the PC, is at the most posterior 
inferior part of V3 and is nearest to the MSP; 

2) checking whether Qi and Q2 have the same location; if yes, V3 is connected 
with V4, and the current connection is terminated; 

3) if no, determining a rectangular area around a pixel on the axial slice next to 
the reference sUce Aj^ in the superior direction, the reference slice being the 
axial slice with Q2; 

4) growing a region on the axial slice next to ^j, in superior direction; if the 
number of pixels grown is zero, go to (7) to force the connection; 

5) checking whether the grown region connects with the subregion 1 of V4; if 
not, connect the grown region wifli flie subregion 1 of V4; 

6) adding the grown region into subregion 1 of V4; assigning the current axial 
slice as the reference slice Aj,; locate the intersection pixel (IP) of the line 
comiecting Qi and Q2 with the axial plane in the superior direction of the 
grown region and make Q2 , the IP, then go to (2); 

7) applying a brute-force connection, by considering the IP in each axial slice as 
the connection. 
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50. An ^paratus arranged to perform a method according to any of the preceding 
claims. 

51. An ^aratus according to claim 50, wherein the ^aratus is a computing 
apparatus. 

52. A computer program product comprismg computer program instructions readable 
by a computer apparatus to cause the computer q)paratus to perform a method 
according to any of the preceding claims. 

53. A method of quantifying the ventricular system by counting the number of voxels 
of the ventricular system extracted by a method as claimed in any preceding claim and 
midtiplying this count by voxel volume. 



